Bounds on Quantum Correlations in Bell Inequality Experiments 
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Bell-inequality violation is one of the most widely known manifestations of entanglement in quan- 
tum mechanics; indicating that experiments on physically separated quantum mechanical systems 
cannot be given a local realistic description. However, despite the importance of Bell inequalities, 
it is not known in general how to determine whether a given entangled state will violate a Bell 
inequality. This is because one can choose to make many different measurements on a quantum sys- 
tem to test any given Bell inequality and the optimization over measurements is a high-dimensional 
variational problem. In order to better understand this problem we present algorithms that pro- 
vide, for a given quantum state and Bell inequality, both a lower bound and an upper bound on 
the maximal violation of the inequality. In many cases these bounds determine measurements that 
would demonstrate violation of the Bell inequality or provide a bound that rules out the possibility 
of a violation. Both bounds apply techniques from convex optimization and the methodology for 
creating upper bounds allows them to be systematically improved. Examples are given to illustrate 
how these algorithms can be used to conclude definitively if some quantum states violate a given 
Bell inequality. 

PACS numbers: 03.65.Ud, 03.67.Mn, 02.60.Pn 
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I. INTRODUCTION 

The nature of entanglement has always been be- 
wildering ever since its first appearance in the litera- 
ture [l|, Q . This is more so after Bell's seminal work 
in 1964 [3], in which he showed, using what is now known 
as a Bell inequality, that some experimental statistics of 
the spin-singlet state are intrinsically incompatible with 
local realism. 

For a long time after that, it seems to have been gen- 
erally assumed that entanglement and the violation of 
Bell inequalities are synonymous. The first counterex- 
ample to that commonly held intuition was provided by 
Werner Q, where he showed that there are bipartite 
mixed two-qudit states, now known as the Werner states, 
that are entangled and yet do not violate a large class of 
Bell inequalities. 

Soon after that, it was demonstrated by Gisin [5| and 
later by Gisin and Peres Q that all bipartite pure en- 
tangled states violate the Bell-Clauser-Horne-Shimony- 
Holt inequality 0] ; the generalization to multipartite pure 
entangled states was also carried out by Popescu and 
Rohrlich by invoking appropriate postselection [f| . Three 
years later, Horodecki et al. provided the first analytic 
criterion [9( to determine if a two-qubit state violates 
the Bell-Clauser-Horne-Shimony-Holt (henceforth abbre- 
viated as Bell-CHSH) inequality. 

The Horodecki criterion is, unfortunately, also the 
only analytic criterion that we have in determining if a 
broad class of quantum states, namely two-qubit states, 
can be simulated classically. For specific quantum states, 
there are examples where explicit local hidden variables 
(LHV) models have been constructed to reproduce the 
quantum mechanical predictions, thereby ruling out the 
possibility that these quantum states may violate a Bell 
inequality 0,[I3|. 
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In general, however, to determine if a quantum state 
violates a Bell inequality is a high-dimensional variational 
problem, which requires a nontrivial optimization of a 
Hermitian operator B (now known as the Bell opera- 
tor (Tl| ) over the various possible measurement settings 
that each observer may perform. This optimization does 
not appear to be convex and is possibly AP-hard 

Except for the simplest scenario where one deals with 
Bell-CHSH inequality [7j, in conjunction with a two-qubit 
state , or a maximally entangled pure state |j, [l3[ , and 
its mixture with the completely mixed state [lj] , very few 
analytic results for the optimal measurements are known. 
As such, for the purpose of characterizing quantum states 
that are incompatible with local realistic description, ef- 
ficient algorithms to perform this state-dependent opti- 
mization are very desirable. 

Bell-inequality violation is also relevant in various 
aspects of quantum information processing, in particu- 
lar, quantum teleportation [l5j |. quantum key distribu- 
tio n fltjl \v?\ , and reduction of communication complex- 
ity [18| ■ Recently, it has even been argued [17| that Bell- 
inequality violation is necessary to guarantee the secu- 
rity of some entanglement-based quantum key distribu- 
tion protocols. 

On the other hand, state-independent bounds of 
quantum correlations have also been investigated since 
the early 1980's. In particular, Tsirelson [n| has demon- 
strated, using what is now known as Tsirelson's vec- 
tor construction, that in a Bell-CHSH setup, bipartite 
quantum systems of arbitrary dimensions cannot exhibit 
correlations stronger than 2^/2 - a value now known as 
Tsirelson's bound. Recently, analogous bounds for more 
complicated Bell inequalities have also been investigated 
by Filip p a nd Svozil pfjj . Buhrman and Massar [2l| . 
Wehner |23|. Toner [23|], Avis et al. [24| and Navascues 
et al. [25]. On a related note, bounds on quantum 
correlations for given local measurements, rather than 
given quantum state, have also been investigated by Ca- 
bello Hi] and Bovino et al. ^2l\ . 

In this paper, we will present, respectively, in Sec. Ill Bl 
and Sec. Ill Cl two algorithms that were developed to pro- 
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vide a lower bound and an upper bound on the max- 
imal expectation value of a Bell operator for a given 
quantum state. The second algorithm is another in- 
stance where a nonlinear optimization problem is ap- 
proximated by a hierarchy of semidefinite programs, each 
giving a better bound of the original optimization prob- 
lem [25|, [2^, [3(|. In its simplest form, it provides 
a bound that is apparently state- independent, and thus 
provides a (not necessarily tight) bound on the maximum 
attainable quantum correlations in a Bell inequality ex- 
periment. 

In Sec. IIHl we will derive, based on the second al- 
gorithm, a necessary condition for a class of two-qudit 
states to violate the Bell-CHSH inequality. We will also 
demonstrate how the two algorithms can be used in tan- 
dem to determine if some quantum states violate a Bell 
inequality. Some limitations of these algorithms will then 
be discussed. We will conclude with a summary of results 
and some possibilities for future research. 

Throughout, boldfaced Latin letters, e.g., x will be 
used to denote a column vector whereas and 1, re- 
spectively, represents a zero block matrix and an identity 
matrix. Moreover, the (i, j) entry of a matrix M will be 
denoted by [M]y. 

II. BOUNDS ON QUANTUM CORRELATIONS 
A. Preliminaries 

Bell inequalities are inequalities derived from the as- 
sumption of a general LHV model. A particular Bell 
inequality deals with a specific experimental setup, say 
a source that distributes pairs of particles to two experi- 
menters (hereafter called Alice and Bob) , and where each 
of them can perform, respectively, vha and tub alterna- 
tive measurements that would each generate ha and ns 
distinct outcomes. For each of these setups, a Bell in- 
equality places a bound on the experimental statistics ob- 
tained from the corresponding Bell experiments. If there 
exists a LHV model that saturates the bound then the 
inequality is said to be tight [32j|. In what follows, we will 
adopt the notation introduced in Ref. [3l| and refer to a 
tight Bell inequality for such an experimental setup as a 
Bell- toa wbi-a^b inequality. 

A Bell inequality for correlations, such as the Bell- 
CHSH inequality pj typically involves statistical con- 
straints on some linear combination of correlation func- 
tions. Similarly, a Bell inequality for probabilities, such 
as the Bell-Clauser-Horne (henceforth abbreviated as 
Bell-CH) inequality [33j|, places bounds on some linear 
combinations of joint and marginal probabilities of exper- 
imental outcomes. In either case, a general Bell inequality 
takes the form: 

Slhv < /3lhv, (1) 

where /3lhv is a re al number and 5lhv involves a specific 
linear combination of correlation functions or joint and 
marginal probabilities of experimental outcomes. 

To compare with predictions given by quantum me- 
chanics, these correlation functions, or probabilities, 
are calculated using the quantum mechanical rules. 



The bounds on 5lhv then translate into corresponding 
bounds /3lhv on the expectation value of some Hermitian 
observable that describes the Bell inequality experiment, 
this observable is known as the Bell operator B [Tl] . The 
restriction that the given Bell inequality is satisfied in the 
experiment is then 

S QM (p,B) = tr(pB)<f3 LHV . (2) 

The Bell operator depends on the choice of measure- 
ments at each of the sites (polarizer angles for exam- 
ple). These measurements will be described by a set of 
Hermitian operators {O m }. For correlation inequalities 
these are simply the measured observables at each stage 
of the Bell measurement, while for general probability in- 
equalities the O m are elements of the positive-operator- 
valued measures (POVMs) that describe the measure- 
ments at each site. We will denote this expectation value 
by Sqm(p, {O m }) when we want to emphasize its depen- 
dence on the choice of local Hermitian observables O m . 
Ideally the choice of measurement should give the maxi- 
mal expectation value of the Bell operator, for which we 
will give the notation 

Sqm(p) = max Sq M (/0, {O m }). (3) 

It is this implicitly-defined function that will give us infor- 
mation about which states violate a given Bell inequality. 

As an example, let us recall the Bell-CHSH inequal- 
ity 0, which is a dichotomic (i.e., two-outcome) Bell cor- 
relation inequality that involves two parties, and where 
the two possible measurement outcomes are assigned the 
values ±1: 

Slhv = E{A 1 ,B 1 ) +E{A 1 ,B 2 ) + E{A 2 ,B 1 ) - E(A 2 , B 2 ) 
< 2- (4) 

In the above expression, the correlation function 
E(Ak,Bi) represents the expectation value of Alice's 
measurement outcome times Bob's measurement out- 
come, given that she has chosen to measure the observ- 
able Ak and he has chosen to measure the observable 
B[. In quantum mechanics, these correlation functions 
are computed using 

E(A k ,B l )=tx(pA k ®B l ). (5) 

Substituting this into Eq. (g]) and comparing with Eq. 
one finds that the corresponding Bell operator reads 

B c ^=A l ®{B l +B 2 ) + A 2 ®{B 1 -B 2 ). (6) 

To determine the maximal Bell-inequality violation 
for a given p, <Sqm (p) , requires a maximization by varying 
over all possible choices of {O m }, i.e., Ak and Bi in the 
case of Eq. ([6]) . Whether we are interested in correlation 
inequalities or in Bell inequalities for probabilities the 
(bipartite) Bell operator has the general structure 

B = J2 b KLA K ®B L . (7) 

K,L 

In the case of a Bell inequality for probabilities the indices 
K, L are collective indices describing both a particular 
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measurement setting and a particular outcome for each 
observer. For correlation inequalities they refer simply 
to the measurement settings as in the Bell-CHSH case 
described in detail above. 

In what follows, we will present two algorithms which 
we have developed specifically to perform the maximiza- 
tion over choice of measurements. The first, which we will 
abbreviate as LB, provides a lower bound on the maxi- 
mal expectation value and can be implemented for any 
Bell inequality. This bound makes use of the fact that 
the objective function <Sqm is bilinear in the observables 
O m , that is it is linear in the Ak for fixed Bl and like- 
wise linear in the Bl for fixed Ak- The second bound, 
which we will abbreviate as UB, provides an upper bound 
on 5qm(p) by regarding Sqm(p, {O m }) as a polynomial 
function of the variables that define the various O m and 
applying general techniques for finding such bounds on 
polynomials HEHi]. 

Both of these make use of convex optimization tech- 
niques in the form of a semidefinite program (SDP) [3~il 
[35j . A semidefinite program is an optimization over Her- 
mitian matrices. The objective function depends linearly 
on the matrix variable (as expectation values do in quan- 
tum mechanics for example) and the optimization is car- 
ried out subjected to the constraint that the matrix vari- 
able is positive semidefinite and satisfies various affine 
constraints. Any semidefinite program may be written in 
the following standard form: 

maximize — tr [i*b^] , (8a) 
subject to tr [F m Z] = c m V m, (8b) 
Z > 0, (8c) 

where Fq and all the F m 's are Hermitian matrices and the 
c m are real numbers that together specify the optimiza- 
tion; Z is the Hermitian matrix variable to be optimized. 

An SDP also arises naturally in the inequality form, 
which seeks to minimize a linear function of the opti- 
mization variables x £ R™, subjected to a linear matrix 
inequality: 

minimize x T c' (9a) 
subject to Go + ^x m G m > 0. (9b) 

m 

As in the standard form, Go and all the G m 's are Hermi- 
tian matrices, while c' is a real vector of length n. 



B. Lower Bound on >Sqm(p) 



Back in 2001, Werner and Wolf [37} presented a simi- 
lar iterative algorithm, by the name of See-Saw iteration, 
to maximize the expectation value of the Bell operator for 
a correlation inequality involving only dichotomic observ- 
ables [39( . As a result we will focus here on the (straight- 
forward) generalization to the widest possible class of 
Bell inequalities. In the work of Werner and Wolf [37[ 
it turned out that once the dichotomic observables for 
one party are fixed, optimization of the other party's ob- 
servables can be carried out explicitly. This turns out 
to be true for any dichotomic Bell inequality and we will 
return to this question in Sec. Ill B 31 

1. General Settings 

Firstly we must develop a more explicit notation for a 
general Bell inequality for probabilities [3a |. Let us con- 
sider a HeW-mAmBriAnB inequality for probabilities. We 
will denote the POVM element associated with the K th 
outcome of Alice's fc th measurement by A% while Bf is 
the POVM element associated with the A th outcome of 
Bob's I th measurement. Moreover, let dA and d B , respec- 
tively, be the dimension of the state space that each of 
the A% and Bf acts on. Then it follows from Born's rule 
that 

P f B {k,l)=tv(pAt®Bt) (10a) 

p1(k) = tr (p Al ® l dB ) , p B (I) = tr (p l dA ® B?) , 

(10b) 

where p\\ (k> re fe rs to the joint probability that the K th 
experimental outcome is observed at Alice's site and the 
A th outcome at Bob's, given that Alice performs the fc th 
and Bob performs the / th measurement. The marginal 
probabilities p^(fe) and p B (l) are similarly defined. A 
general Bell operator for probabilities can then be ex- 
pressed as 

m A n A m B n B 

k=l K=l 1 = 1 A=l 

where 6£ ; A are determined from the given Bell inequal- 
ity. Note that the sets of POVM elements {Ag}"^ and 
{B?}™! satisfy 

"A n B 

Y t A% = U A and J2 B i= ld B V k,l, (12a) 

K=l A=l 

Afe>0, B?>0 V k,l, K ,X. (12b) 



The key idea behind the LB algorithm is to realize 
that when measurements for all but one party are fixed, 
the optimal measurements for the remaining party can be 
obtained efficiently using convex optimization techniques, 
in particular an SDP. Thus we can fix Bob's measure- 
ments and find Alice's optimal choice, at least numer- 
ically; with this optimized measurements for Alice, we 
can further find the optimal measurements for Bob (for 
this choice of Alice's settings), and then Alice again and 
so on and so forth until iSqm (p,{Om}) converges within 
the desired numerical precision [36j . 



2. Iterative Semidefinite Programming Algorithm 

To see how to develop a lower bound on Sqm(p) 
by fixing the observables at one site and optimizing the 
other, we observe that upon substituting Eq. (jTT|) into 
Eq. ((2|) , the left-hand-side of the inequality can be rewrit- 
ten as 

S QM (p, A K k ,Bf) = ]T tr (p B , Bf) , (13) 
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where p B x = J2k, K b M tr A [p {A% ® U B )}- ■ is the 
partial trace over subsystem A. 

Notice that if p B f are held constant by fixing all of 
Alice's measurement settings (given by the set of A?) 
then ps ; is a constant matrix independent of the Bf. 
Thus the objective function is linear in these variables. 
The constraints that {-B; A }"2i form a POVM for each 
value of I is a combination of affine and matrix nonneg- 
ativity constraints. As a result it is fairly clear that the 
following problem is an SDP in standard form 



maximize{ B \ j <Sqm (p, A% , B A ) 

»B 

subject to 2_. = ld B V I, 



A=l 



B ; A >0 V l,X 



(14a) 
(14b) 
(14c) 



Explicit forms for the matrices F m and values c m of 
Eq. ([8]) can be found in Appendix |A"1 

Exactly the same analysis follows if we fix Bob's 
measurement settings, and optimize over Alice's POVM 
elements instead. To arrive at a local maximum of 
Sqm(p, A%, Bf), it therefore suffices to start with some 
random measurement settings for Alice (or Bob) , and op- 
timize over the two parties' settings iteratively. A (non- 
trivial) lower bound on 5qm(p) can then be obtained by 
optimizing the measurement settings starting from a set 
of randomly generated initial guesses. 

It is worth noting that in any implementation of 
this algorithm, physical observables {A%,B^} achieving 
the lower bound are constructed when the correspond- 
ing SDP is solved. In the event that the lower bound 
is greater than the classical threshold /3lhv, then these 
observables can, in principle, be measured in the labo- 
ratory to demonstrate a Bell-incquality violation of the 
given quantum state. 

Wc have implemented this algorithm in MAT- 
LAB to search for a lower bound on Sqm(p) in the 
case of Bell-CH [33|, 13322, ^4422, ^2233 and J 2 244 inequal- 
ity [3ll | . and with the local dimension d up till 32. Typ- 
ically, with no more than 50 iterations, the algorithm 
already converges to a point that is different from a local 
maximum by no more than 10~ 9 . To test against the 
effectiveness of finding Sqm(p) using LB, we have ran- 
domly generated 200 Bell-CH violating two-qubit states 
and found that on average, it takes about 6 random 
initial guesses before the algorithm gives <Sqm(p, {O m }) 
that is close to the actual maximum (computed using 
Horodecki's criterion within 10 -5 . Specific examples 
regarding the implementation of this algorithm will be 
discussed in Sec. IIIII 

Two other remarks concerning this algorithm should 
now be made. Firstly, the algorithm is readily general- 
ized to multipartite Bell inequalities for probabilities: one 
again starts with some random measurement settings for 
all but one party, and optimizes over each party itera- 
tively. Also, it is worth noting that this algorithm is not 
only useful as a numerical tool, but for specific cases, 
it can also provide useful analytic criterion. In particu- 
lar, when applied to the Bell-CH inequality [33| for two- 
qubit states, the analysis for dichotomic observables dis- 
cussed in the next subsection allows one to recover the 
Horodecki's criterion [9[, i.e., the necessary and sufficient 



condition for two-qubit states to violate the Bell-CHSH 
inequality 0,[4l|. 



3. Two-outcome Bell Experiment 

We will show that, just as in the case of correlation 
inequalities [37| , the local optimization can be solved ana- 
lytically for two-outcome measurements. If we denote by 
"±" the two outcomes of the experiments, it follows from 
Eq. (|12p that the POVM element Bf can be expressed 
as a function of the complementary POVM element B^~ , 



i.e., Bf = 1 
then have 



B+, subjected to < B+ < l a 



We 



A=± 



tr(p B xS ; A ) =tr[(p B +-p Br )B] 



tr 



Pb; 



The above expression can be maximized by setting the 
positive semidefinite operator B~f to be the projector onto 
the positive eigenspace of p B + —p B -. In a similar manner, 
we can also write 



tr (p B > B? 

x=± 



= tr 



Pbt ~ Pb+ 



b; 



tr 



Pb+ 



which can be maximized by setting B t to be the projec- 
tor onto the nonpositive eigenspace of p B + — p B - . Notice 

that this choice is consistent with our earlier choice of B f 
for the + outcome POVM element in that they forms a 
valid POVM. Since there can be no difference in these 
maxima, we may write the maximum as their average, 
i.e., 



A=± 



tr 



Pl;> ) = ~ 



Pb+ ~ Pbt 



1 



tr 



Pb? 



where ||0|| is the trace norm of the Hcrmitian operator 

[42j. Carrying out the optimization for each of the 

1 settings, the optimized 5qm(p, A%, Si ), as an implicit 
function of Alice's POVM {A%}, is given by 



It. 



\Pb+ - Pb; 



l A=± 



Pb?- 



(15) 



An immediate corollary of the above result is that 
for the optimization of a two-outcome Bell operator for 
probabilities, it is unnecessary for any of the two ob- 
servers to perform generalized measurements described 
by a POVM; von Neumann projective measurements 
are sufficient (43|. In practice, this simplifies any ana- 
lytic treatment of the optimization problem as a generic 
parametrization of a POVM is a lot more difficult to deal 
with, thereby supporting the simplification adopted in 
Ref. [U. 

Nevertheless, it may still be advantageous to consider 
generic POVMs as our initial measurement settings when 
implementing the algorithm numerically. This is because 
the local maximum of Sqm(p, {O m }) obtained using the 
iterative procedure is a function of the initial guess. In 
particular, it was found that the set of local maxima at- 
tainable could change significantly if the ranks of the ini- 
tial measurement projectors are altered. As such, it seems 
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necessary to step through various ranks of the starting 
projectors to obtain a good lower bound on Sqm (p) ■ Even 
then, we have also found examples where this does not 
give a lower bound on Sqm(p) that is as good as when 
generic POVMs are used as the initial measurement op- 
erators. 



C. Upper Bound on Sqm(p) 

A major drawback of the above algorithm, or the 
analogous algorithm developed by Werner and Wolf [13] 
for Bell correlation inequalities is that, except in some 
special cases, it is generally impossible to tell if the max- 
imal Sqm(p, {O m }) obtained through this optimization 
procedure corresponds to the global maximum Sqm (p) ■ 

Nontrivial upper bounds on iSqm(p), nevertheless, 
can be obtained by considering relaxations of the global 
optimization problem. In a relaxation, a (possibly non- 
convex) maximization problem is modified in some way 
so as to yield a more tractable optimization that bounds 
the optimization of interest. One example of a variational 
upper bound that exists for any optimization problem is 
the Lagrange dual optimization that arises in the method 
of Lagrange multipliers (35[ . 

To see how to apply existing studies in the optimiza- 
tion literature to find upper bounds on Sqm(p), let us first 
remark that the global objective function Sqm{p, {O m }) 
can be mapped to a polynomial function in real variables, 
for instance, by expanding all the local observables {O m } 
and the density matrix p in terms of Hermitian basis oper- 
ators. In the same manner, matrix equality constraints, 
such as that given in Eq. (|12ap can also mapped to a 
set of polynomial equalities by requiring that the matrix 
equality holds component wise. Now, it is known that a 
hierarchy of global bounds of a polynomial function, sub- 
jected to polynomial equalities and inequalities, can be 
achieved by solving suitable SDPs [28|, [2^| . Essentially, 
this is achieved by approximating the original nonconvex 
optimization problem by a series of convex ones in the 
form of a SDP, each giving a better bound of the original 
polynomial objective function. 

At the bottom of this hierarchy is the lowest order 
relaxation provided by the Lagrange dual of the original 
nonconvex problem. By considering Lagrange multipli- 
ers that depend on the original optimization variables, 
higher order relaxations to the original problem can be 
constructed to give tighter upper bounds on Sqm (p) (see 
Appendix FBI for more details). 

In the following, we will focus our discussion on a 
general two-outcome Bell (correlation) inequality, where 
the observables {O m } are only subjected to matrix equal- 
ities. In particular, we will show that the global optimiza- 
tion problem for these Bell inequalities is a quadratically- 
constrained quadratic-program (QCQP), i.e., one whereby 
the objective function and the constraints are both 
quadratic in the optimization variables. Then, we will 
demonstrate explicitly how the Lagrange dual of this 
QCQP, which is known to be an SDP, can be con- 
structed. The analogous analytic treatment is apparently 
formidable for higher order relaxations. Nonetheless, 
there exists third-party MATLAB [4(| toolbox known as 
the SOSTOOLS which is tailored specifically for this kind 



of optimization problem [ifl |46| . 

Numerically, we have implemented the algorithm for 
several two-outcome correlation inequalities and will dis- 
cuss the results in greater detail in Sec. IIII1 For a general 
Bell inequality where each O m is also subjected to a lin- 
ear matrix inequality (LMI) like Eq. (| 1 2b[) , the algorithm 
can still be implemented, for instance, by requiring that 
all the principle submatrices of O m have nonnegative de- 
terminants 42]. This then translates into a set of poly- 
nomial inequalities which fit into the framework of a gen- 
eral polynomial optimization problem (see Appendix IB"]) . 
However a more effective approach would retain the struc- 
ture of linear matrix inequalities constraining a polyno- 
mial optimization problem; we leave the investigation of 
these bounds to further work. 



1. Global Optimization Problem 

Now, let us consider a dichotomic Bell correlation 
inequality where Alice an Bob can respectively perform 
rriA and ms alternative measurements. A general Bell 
correlation operator for such an experimental setup can 
be written as [47| 

m.A mB 
k=\ l=\ 

where b^i are determined from the given Bell correla- 
tion inequality, Ok for k = 1, . . . , refers to the k th 
Hermitian observable measured by Alice, and Oi+ mA for 
I = 1, . . . , vriB refers to the ^ th Hermitian observable mea- 
sured by Bob. Furthermore, by convention, these di- 
chotomic observables can be chosen to have eigenvalues 
±1 [H| and thus 

OlO m = (O m ) 2 = l d V m, (17) 

where we have assumed for simplicity that all the local 
observables O m act on a state space of dimension d [49| . 

The global optimization problem derived from a di- 
chotomic Bell correlation inequality thus takes the form 
of 

maximize tr (p B) , (18a) 
subject to 2 m = l d V m — 1,2, .... (18b) 

For any m x n complex matrices, we will now define 
vec(A) to be the m ■ n dimensional vector obtained by 
stacking all columns of A on top of one another. By col- 
lecting all the vectorized observables together 

w f = [vec(Oi) f vec(0 2 ) f •■■ vec{O mA+mB ^}, 

and using the identity 

tr(p O k ® O l+mA ) = vec{O k y(Vp) TA vec(O l+mA ), (19) 

with V being the flip operator such that V^ji) j j) = \j)\i) 
and ( ) Ta being the partial transposition with respect to 
subsystem A, we can write the objective function more 
explicitly as 

Sq M (p, {O m }) = tr(p B) = -wtft w (20) 
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where 

n _ 1 ( -b®R\ 

b is a rriA x tub matrix with [b]ki — bki [c.f. Eg.(|16[)] and 
R = (V p) Ta ■ In this form, it is explicit that the objective 
function is quadratic in vec(O m ). Similarly, by requiring 
that the matrix equality holds component-wise, we can 
get a set of equality constraints, which are each quadratic 
in vec(O m ). The global optimization problem (fl8|) is thus 
an instance of a QCQP. 

On a related note, for any Bell inequality experiments 
where measurements are restricted to the projective type, 
the global optimization problem is also a QCQP. To see 
this, we first note that the global objective function for 
the general follows from Eq. © and Eq. fl7J), 

is always quadratic in the local Hermitian observables 
{Ak, Bl}- The requirement that these measurement op- 
erators are projectors amounts to requiring 

A 2 K = A K , B\=B U V K, L, (21) 

which are quadratic constraints on the local Hermitian 
observables. Since we have shown in Sec. Ill B 31 that 
for any two-outcome Bell inequality for probabilities, it 
suffices to consider projective measurements in optimiz- 
ing Sqm(p, {O m }), it follows that the global optimization 
problem for these Bell inequalities is always a QCQP. 

2. State-independent Bound 

As mentioned above, the lowest order relaxation to 
the global optimization problem (|18[) is simply the La- 
grange dual of the original QCQP. This can be obtained 
via the Lagrangian of the global optimization problem 

m 

mA+rriB 

C({O m }, A ro ) = 5 qm (p, {O m })~ tr [ A ™ (°m - 

771—1 

where A m is a matrix of Lagrange multipliers associated 
with the m th matrix equality constraint. With no loss of 
generality, we can assume that A m 's are Hermitian. 

Notice that for all values of {O m } that satisfy 
the constraints, the Lagrangian C(p, {O m },A m ) = 
Sqm(p, {Om}) ■ As a result, if we maximize the La- 
grangian without regard to the constraints we obtain an 
upper bound on the maximal expectation value of the 
Bell operator 

max£(p, {O m }, A m ) > 5 QM (/o). (23) 

The Lagrange dual optimization simply looks for the best 
such upper bound. 

In order to maximize the Lagrangian we rewrite the 
Lagrangian with Eq. ([20]) and the identity 

tr (A rn O m Ol n ) = vec(O m ) t ld <8> A m vec(O m ), (24) 

to obtain 

mA+rriB 

£(w,A m ) = -wfj]w+ trA ™> ( 25 ) 

m— 1 



where Q = Q Q + 0™=| ms U ® A m . Note that each of 
the diagonal blocks Id® A m is of the same size as the 
matrix R. 

To obtain the dual optimization problem, we maxi- 
mize the Lagrangian over w to obtain the Lagrange dual 
function 

g(A m ) = sup£(w, A m ). (26) 

w 

As noted above g(A m ) > Sqm(p) for all choices of A m . 
Moreover, this supremum over w is unbounded above un- 
less ft > 0, in which case the supremum is attained by 
setting w = in Eq. (|25[) . Hence, the Lagrange dual 
optimization, which seeks for the best upper bound of 
Eq. (|18p by minimizing Eq (|26[) over the Lagrange multi- 
pliers, reads 

mA+m B 

minimize tr A m , 

4=1 

subject to fl > 0. (27) 

By expanding A m in terms of Hermitian basis operators 
satisfying Eq. (|A1[) . 

d 2 -i 

A m = y A mn cr n , (28) 

n=Q 

the optimization problem (|27p is readily seen to be an 
SDP in the inequality form ©. 

For Bell-CHSH inequality and the correlation equiv- 
alent of /3322 inequality, it was observed numerically that 
the upper bound obtained via the SDP (|27p is always 
state-independent. In fact, for 1000 randomly gener- 
ated two-qubit states, and 1000 randomly generated two- 
qutrit states, the Bell-CHSH upper bound of Sqm(p) ob- 
tained through (|27D was never found to differ from the 
Tsirelson bound [lj| by more than 10 -7 . In fact by find- 
, ing an explicit feasible solution to the optimization prob- 
lem dual to Eq. (|2T|) . Wehner has shown that the upper 
bound obtained here can be no better than that obtained 
by Tsirelson's vector construction for correlation inequal- 
ities [50l |. 

In a similar manner, we have also investigated the 
upper bound of Sqm(p) for some dichotomic Bell prob- 
ability inequalities using the lowest order relaxation to 
the corresponding global optimization problem. Interest- 
ingly, the numerical upper bounds obtained from the ana- 
log of Eq. (|27p for these inequalities, namely the Bell-CH 
inequality, the 13322 inequality and the -Z4422 inequality, 
are also found to be state-independent and are given by 
0.2071067, 0.375 and 0.6693461 respectively. 

3. State- dependent Bound 

Although the state-independent upper bounds ob- 
tained above are interesting in their own right, our main 
interest here is to find an upper bound on 5qm(p) that 
does depend on the given quantum state p. This can be 
obtained, with not much extra cost, from the Lagrange 
dual to a more-refined version of the original optimization 
problem. 
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To appreciate that, let us first recall that each di- 
chotomic Hermitian observable O m can only have eigen- 
values ±1. It follows that their trace 



%rn — tr(O m ), 
can only take on the following values 



= -d,-d + 2, 



,d-2,d. 



(29a) 



(29b) 



In particular, if z m — ±d for any m, then O m — ±l<j and 
it is not difficult to show that the Bell-CHSH inequality 
cannot be violated for this choice of observable. 

Better Lagrange dual bounds arise from taking these 
additional constraints (|29|) explicitly into account. We 
found it most convenient to express the original opti- 
mization problem in terms of real variables given by the 
expansion coefficients of O m in terms of a basis for Her- 
mitian matrices that includes the (traceless) Gcll-Mann 
matrices and the identity matrix [c.f. Eq. (jAlj) ]. For de- 
tails see Appendix O The result is a set of SDPs, one 
for each of the various choices of z m . The lowest order 
upper bound on Sqm(p) can then be obtained by step- 
ping through the various choices of z m given in Eq. (|29b|) , 
solving each of the corresponding SDP, and taking their 
maximum. The results of this approach will be discussed 
later, for now it suffices to note that tighter bounds can 
be obtained which are explicitly state dependent. 



4- Higher Order Relaxations 

The higher order relaxations simply arise from allow- 
ing the Lagrange multipliers A to be polynomial func- 
tions of {O m } rather than constants. In this case, it is 
no longer possible to optimize over the primal variables 
in the Lagrangian analytically but let us consider the fol- 
lowing optimization 



minimize 



7 



subject to 7 - Sqm(p, x) = /u(x) + Ai(x)/ eq ,i(x), 

(30) 

where each of the A^ 's are polynomial functions of x and 
/z(x) is a sum of squares polynomial (SOS) and there- 
fore nonnegative. That is fi(x) — J2j[hj( x )] 2 > for 
some set of real polynomials hj(x). The variables of 
the optimizations are x and the coefficients that de- 
fine the polynomials /x(x) and A^. Notice that we have 
7 > Sqm (p, x) whenever the constraints arc satisfied 
so that once again we have a global upper bound on 
<Sqm (p, x) . This optimization can be implemented nu- 
merically by restricting /i(x) and A^ to be of some fixed 
degree. The Lagrange dual optimization (|27[) arises from 
choosing the degree of Xi to be zero. It is known that 
for any fixed degree this optimization is an SDP [28|, [2i| 
and we have implemented up to degree four using SOS- 
TOOLS [45|,[46|. Schmiidgen's theorem guarantees that 
by increasing the degree of the polynomials in the relax- 
ation we obtain bounds approaching the true maximum 
Sqm(p)- This is a special case of the general procedure 
described in [28|, |4y, |5I( which is able to handle inequality 
constraints. For more details see Appendix iBl 



III. APPLICATIONS & LIMITATIONS 

In this section, we will look at some concrete exam- 
ples of how the two algorithms can be used to determine 
if some quantum states violate a Bell inequality. Specif- 
ically, we begin by looking at how the second algorithm 
can be used to determine, both numerically and analyti- 
cally, if some bipartite qudit state violates the Bell-CHSH 
inequality. Then in Sec. IIIIB1 we demonstrate how the 
two algorithms can be used in tandem to determine if 
a class of two-qubit states violate a Bell-3322 inequal- 
ity [3l[. We will conclude this section by pointing out 
some limitations of the UB algorithm that we have ob- 
served. 



A. Bell-CHSH violation for Two-Qudit States 

The Bell-CHSH inequality, as given by Eq. is one 
that amounts to choosing 



I I 

I -1 



(31) 



For low-dimensional quantum systems, an upper bound 
on Sqm(p) can be efficiently computed in MATLAB fol- 
lowing the procedures described in Sec. Ill C 31 However, 
for high-dimensional quantum systems, intensive com- 
putational resources are required to compute this upper 
bound, which may render the computation infeasible in 
practice. Fortunately, for a specific class of two-qudit 
states, namely those whose coherence vectors [52| van- 
ish, it can be shown that (Appendix |D| their 5qm(p) = 
maxg CHSH (Schsh)p cannot exceed 



max 



2v / 2sid 



\ 



n 



2d 2 - 



■2i-l 



"2i 



2d 2 



k,l 



,( 32 ) 

where si is the largest singular value of the matrix R' 
defined in Eq. (|C2a|) . and z m is the trace of the dichotomic 
observable O m given in Eq. (|2U| . 

To violate the Bell-CHSH inequality, we must have 
Sqm(p) > 2, hence for this class of quantum states, the 
Bell-CHSH inequality cannot be violated if 



max 



2d 2 



2d 2 



k.i 



ZkZl+2 , 1 

< 1. 



(33) 

Since the Bell-CHSH inequality is tight Eq. (f33|) guar- 
antees the existence of a LHV model for this experi- 
mental setup [Eij]. Essentially, the semianalytic upper 
bound Eq. (|32[) was obtained by considering a particu- 
lar choice of Lagrange multipliers in the Lagrange dual 
function ((26|) . Hence, it is generally not as tight as the 
upper bound obtained numerically using the procedures 
described in Sec. Ill C 31 (see Appendix iDl for details). 

As an example, consider the d-dimensional isotropic 
state, i.e., a mixture of the d-dimensional maximally en- 
tangled state | VP J) = l/v^X^=i an( l tne com- 
pletely mixed state: 



pi* =p|tt£X*3l + (i-p) 



id 2 

d 2 ■ 



(34) 
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As can be verified using the Positive Partial Transposi- 
tion (PPT) criterion [EBI, [EH], this state is entangled if 
and only if p > pEnt = l/(d + 1). Using the procedures 
outlined in Sec. Ill C 31 we can numerically compute, up 
till d = 5, the threshold value of p below which there can 
be no violation of the Bell-CHSH inequality; these critical 
values, denoted by puB-numericai can be found in column 4 
of TableJU Similarly, we can numerically compute the cor- 
responding threshold values, denoted by puB-scmianaiytic, 
using Eq. ([33)1 . It is worth noting that these threshold 
values, as can be seen from column 3 and 4 of Table HI 
agree exceptionally well, thereby suggesting that the com- 
putationally feasible criterion given by Eq. (|33[) may be 
exact for the isotropic states. 



TABLE I: The various threshold values for isotropic states. 
The first column of the table is the dimension of the local 
subsystem d. From the second column to the fifth column, 
we have respectively the value of p beyond which the state is 
entangled pEnt, the value of p below which Eq. (|33p is satisfied 
PuB-semianaiytic, the value of p below which the upper bound 
obtained from lowest order relaxation is compatible with Bell- 
CHSH inequality and the value of p beyond which a Bell- 
CHSH violation has been observed using the LB algorithm. 



obtain the corresponding Bell operator for probabilities: 



d 


PEnt 


PUB- semi analytic 


PUB- numerical 


J>LB 


2 


0.33333 


0.70711 


0.70711 


0.70711 


3 


0.25000 


0.70711 


0.70711 


0.76297 


4 


0.20000 


0.65465 


0.65465 


0.70711 


5 


0.16667 


0.63246 


0.63246 


0.74340 


10 


0.09091 


0.51450 




0.70711 


25 


0.03846 


0.36490 




0.71516 


50 


0.01961 


0.26963 




0.70711 



B. /3322-violation for Two-Qubit States 

Next, we look at how the two algorithms can be used 
in tandem to determine if some two-qubit states violates 
a Bell-3322 inequality, in particular the 13322 inequal- 
ity [3l[. This Bell inequality is interesting in that there 
are quantum states that violate this new inequality but 
not the Bell-CH nor the Bell-CHSH inequality. The ana- 
logue of Horodecki's criterion for this inequality is thus 
very desired for the characterization of quantum states 
that are inconsistent with local realistic description. 

To the best of our knowledge, such an analytic crite- 
rion is yet to be found. However, by combining the two 
algorithms presented above, we can often offer a defini- 
tive, yet nontrivial, conclusion about the compatibility 
of a quantum state description with that given by LHV 
models. To begin with, we recall that the ^3322 inequality 
reads: 

5 LH v=pli(l,l)+pl5(l>2)+pii(l,3)+pii(2 > l) 
+p\ B (2, 2) - p+ B (2, 3) + p+ B (3, 1) - p+ B (3, 2) 
-p+(l)-2p B (l)-p B (2)<0. (35) 

Together with Eq. © and Eq. ([IT))) , one can then 



B, 



x',22 



-M 



At 



-A-2 <8 {Bi + B2) — Aj 



B 2 - A 3 <g> B 1 



(36) 



where we have also used Eq. ()12ap to arrive at this form. 
For convenience, we will adopt the notation that 0„ = 
for m = 1,2,3 and = B m _ 3 for m = 4, 5, 6. In 
these notations, the global optimization problem for this 
Bell inequality can be written as 



maximize tr(p #3322) 



subject to {O ri 



Or, 



(37a) 
(37b) 



which is again a QCQP [57|. The lowest order relaxation 
to this problem can thus be obtained by following similar 
procedures as that described in Sec. Ill CI 

To obtain a state- dependent upper bound on <Sqm(/?) 
for this inequality, we have to impose the analogue of 
Eq. ((29b) . i.e., 



tr(Ojj 



0,1,. 



(38) 



for each of the POVM elements. For small d, numerical 
upper bounds on <Sqm(p) can then be solved for using 
SOSTOOLS. As an example, let's now look at how this 
upper bound, together with the LB algorithm, has en- 
abled us to determine if a class of mixed two-qubit states 
violates the ^3322 inequality. 
The mixed two-qubit state 



PCG = p|*2:l)<*2:l 1 + 0- ~ p) 



<P < 1, 
(39) 



can be understood as a mixture of the pure product state 
|0)|1) and the nonmaximally entangled two-qubit state 
|* 2:1 ) = 4= (2|0)|0) + |1)|1)). As can be easily verified 

I L 

using the PPT criterion [55j, this state is entangled for 
< p < 1. In particular, the mixture with p = 0.85 
was first presented in Ref. (3l| as an example of a two- 
qubit state that violates the 13322 inequality but not the 
Bell-CH nor the Bell-CHSH inequality. 







B 



06 0.16 



C 







D 



81 



1 



FIG. 1: Domains of p where the compatibility between local 
realism and quantum mechanical description given by pcG 
were studied via the 73322 inequality. From right to left are 
respectively the domain of p whereby pea is: (D) found to 
violate the 73322 inequality; (C) found to give a lowest order 
upper bound that is compatible with the 73322 inequality; (B) 
found to give a higher order upper bound that is compatible 
with the 73322 inequality; (A) not known if it violates the 73322 
inequality. 

Given the above observation, a natural question that 
one can ask is, at what values of p does pec violate the 
-^3322 inequality? Using the LB algorithm, we have found 
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that for p > 0.83782 (domain D in Fig. [J), p CG vio- 
lates the -Z3322 inequality. As we have pointed out in 
Sec. Ill BL observables that lead to the observed level of 
/3322-violation can be readily read off from the output of 
the SDP. 

On the other hand, through the UB algorithm, we 
have also found that, with the lowest order relaxation, 
the states do not violate this 3-setting inequality for 
0.16023 < p < 0.83625 (domain C in Fig. [[]); with a 
higher order relaxation, this range expands to 0.06291 < 
p < 0.83782 (domain B in Fig. [1]); the next order re- 
laxation is, unfortunately, beyond what the software can 
handle. 

The algorithms alone therefore leave a tiny gap at 
< P 0.06291 (domain A in Fig. [1]) where we could 
not conclude if pec violates the 13322 inequality. Nev- 
ertheless, if we recall that the set of quantum states 
not violating a given Bell inequality is convex and that 
Pcg{p = 0), being a pure product state, cannot violate 
any Bell inequality, we can immediately conclude that 
Pcg with < p < 0.83782 cannot violate the ^3322 in- 
equality As such, together with convexity arguments, 
the two algorithms allow us to fully characterize the state 
Pcg compatible with LHV theories, when each observer 
is only allowed to perform three different dichotomic mea- 
surements. 



0.25 r- 
0.2 [ 
0.15 \ 



I -0.05 \ 

a> . - % ........ - 

& -0.1- \ 
z> \ 

-0.15- v > : - 

-0.2- '**': — : - - 

-0.25 1 — ■ — 1 — ■ — 1 — ■ — 1 — ■ — 1 — ■ — 1 — ■ — 1 — ■ — 1 — ■ — 1 — ■ — 1 — — 
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p 

FIG. 2: (Color online) Numerical upper bound on maximal 
{Bcu) PH obtained from the UB algorithm using lowest order 
relaxation and Eq. (|38[) . The dotted horizontal line is the 
threshold above which no LHV description is possible. 

the pure product state, ph(p = 0) = |^ / p )(4'p| is actu- 
ally the maximal achievable Bcll-CH violation given by a 
quantum state fl9l ]. 



C. Limitations of the UB algorithm 



IV. CONCLUSION 



As can be seen in the above examples, the UB algo- 
rithm does not always provide a very good upper bound 
for <Sqm(p)- In fact, it has been observed that for pure 
product states, the algorithm with lowest order relaxation 
always returns a state-independent bound (the Tsirelson 
bound in the case of Bell-CHSH inequality). As such, for 
mixed states that can be decomposed as a high-weight 
mixture of pure product state and some other entangled 
state, the upper bound given by UB is typically bad. To 
illustrate this, let us consider the following 1-parameter 
PPT bound entangled state [HSl]: 



PH 



"PEnt 



8p+V 8 P +1 



\%){%\, 0<P<1 (40a) 



where 



|* p ) = |2) Ji^lO) + \/^|2) . (40b) 



It can be shown, for example, using the range crite- 
rion [5§| that this state is entangled when < p < 1. It is, 
however, well-known that a bipartite PPT state cannot 
violate the Bell-CH nor the Bell-CHSH inequality (6(J. 

When tested with the UB algorithm using the lowest 
order relaxation, it turned out that some of these upper 
bounds are actually above the threshold of Bell-CH vio- 
lation (see Fig. [2]). In fact, the upper bound obtained for 



In this paper, we have presented two algorithms 
which can be used to determine, respectively, a lower 
bound (LB) and an upper bound (UB) on the maximal 
expectation value of a Bell operator for a given quantum 
state, i.e., <Sqm(/?)- In particular, we have demonstrated 
how one can make use of the upper bound to derive a 
necessary condition for two-qudit states with vanishing 
coherence vectors to violate the Bell-CHSH inequality 

For low dimensional quantum systems, we have also 
demonstrated how one can make use of the two algo- 
rithms to determine, numerically, if the quantum mechan- 
ical prediction is compatible with local realistic descrip- 
tion. On a separate note, these algorithms have also been 
applied to the search of maximal-Bell-inequality- violation 
in the context of collective measurements without posts- 
election [oil ]. 

As with many other numerical optimization algo- 
rithms, the algorithm to determine a lower bound (LB) 
on <Sqm(p), can only guarantee the convergence to a lo- 
cal maximum. The UB algorithm, on the other hand, 
provides an (often loose) upper bound on Sqm{p)- In 
the event that these bounds agree (up till reasonable nu- 
merical precision) , we know that optimization of the cor- 
responding Bell operator using LB has been achieved. 
This ideal scenario, however, is not as common as we 
would like it to be. In particular, the UB algorithm with 
lowest order relaxation has been observed to give rather 
bad bounds for states with a high-weight mixture of pure 
product states (although it appears that we can often 
rule out the possibility of a violation in this situation 
by convexity arguments as in Sec. IIII B|) . A possibility 
to improve these bounds, as suggested by the work of 
Nie et al. |62| , is to incorporate the Karush-Kuhn- Tucker 
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optimality condition as an additional constraint to the 
problem. We have done some preliminary studies on this 
but have not so far found any improvement in the bounds 
obtained but this deserves further study. 

As of now, we have only implemented the UB al- 
gorithm to determine upper bounds on Sqm(p) for di- 
chotomic Bell inequalities. For Bell inequalities with 
more outcomes, the local Hermitian observables are gen- 
erally also subjected to constraints in the form of a LMI. 
Although the UB algorithm can still be implemented for 
these Bell inequalities by first mapping the LMI to a se- 
ries of polynomial inequalities, this approach seems bla- 
tantly inefficient. Future work to remedy this difficulty 
is certainly desirable. 

Finally, despite the numerical and analytic evidence 
at hand, it is still unclear why the lowest order relax- 
ation to the global optimization problem, as described 
in Sec. Ill C 21 seems always gives rise to a bound that is 
state-independent and how generally this is true. Some 
further investigation on this may be useful, particularly 
to determine whether the lowest order relaxation is al- 
ways state independent even for inequalities that are not 
correlation inequalities. If so this could complement the 
methods of Refs. @,[24|,[25| for finding state- independent 
bounds on Bell inequalities. 
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APPENDIX A: EXPLICIT FORMS OF 
SEMIDEFINITE PROGRAMS 

Here, we provide an explicit form for the matrices 
F m and constants c m that define the semidcfinite pro- 
gram (TrjJ. By setting 

Z = B\ © . . . © Bl B © B\ © . . . e Blf B , 

f = - (p B I © .. . © pbT © p B \ © . e p B ^ B ) , 

in Eq. ([5]) , we see that the inequality constraint (15c)) of the 
SDP entails the positive semidefiniteness of the POVM 
elements {-B^},., and hence Eq. (|L4c|) . On the other 
hand, the equality constraint {8b]), together with appro- 
priate choice of F m and c m , ensures that the normaliza- 
tion condition (|14bp is satisfied. 

In particular, each F m is formed from a direct sum of 
Hermitian basis operators. A convenient choice of such 
basis operators is given by the traceless Gell-Mann matri- 
ces, denoted by {cij^J^ 1 , supplemented by ao — ld/Vd 
such that 

tr (a n a n i) = 5 nn ' and tr (a n ) = \ZdS n0 , (Al) 

where d = d B is the dimension of the state space that each 
acts on. A typical F m then consists of ub diagonal 
blocks of <jj at positions corresponding to the ng POVM 



elements {-B; }a=i in Z for a fixed I. For instance, the 
set of F m 



7TI— 1 
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1 < m < d 2 , 



together with c m = yd 6 m i entails the normalization of 

{-^l }a=h Le -> Ea=i b i = Ub] the remaining {m B ~l)d 2 
F m are defined similarly and can be obtained by shifting 
the nonzero diagonal blocks diagonally downward by ap- 
propriate multiples of n B blocks. The SDP thus consists 
of solving Eq. (JSJ for a m B n B d x mgngd Hermitian ma- 
trix Z subjected to d 2 rriB affine constraints. 



APPENDIX B: SEMIDEFINITE RELAXATION 
TO THE GLOBAL OPTIMIZATION PROBLEM 

The global optimization problem, either in the form 
of Eq. (1T51) for a two-outcome Bell correlation inequality, 
or Eq. (|37[) for a two-outcome Bell inequality for proba- 
bilities, is a QCQP. As was demonstrated in Sec. Ill CI an 
upper bound on Sqm (p) can then be obtained by consid- 
ering the corresponding Lagrange Dual. 

More generally, the global optimization problem can 
be mapped to a real polynomial optimization problem: 

maximize / Q bj(y), (Bla) 
subject to /eq,»(y) = 0, i = 1,2, . . ., N eq , (Bib) 
/in 0qj (y) > 0, j = 1, 2, ... , N meq , (Blc) 

where y is a vector of real variables formed by the ex- 
pansion coefficients of local observables {O m } in terms of 
Hermitian basis operators. 

Results from semialgebraic geometry dictate that 
an upper bound for / bj(y) can be computed using 
Positivstellensatz-bascd relaxations (see, for example, 
Ref. [28| and references therein). In particular, 7 will be 
an upper bound on the constrained optimization problem 
(|B1[) if there exists a set of sum of squares (SOS) /ij(y)'s 
(i.e., nonnegative, real polynomials that can be written as 
J2j[hj(y)] 2 w ith hj(y) being some real polynomials of y), 
and a set of real polynomials Vj{y) such that [28, 46H5l| 

7 - Jobj(y) =Mo(y) + ^(y)/eq,j(y) + Mi(y)/iBeq,i(y) 

3 i 

+ Yj Mii,i 2 (y)/ineq,-ii(y)/i I1 eq,-i 2 (y) + • ■ ■ ■ 

(B2) 

The relaxed optimization problem then consists of min- 
imizing 7 subjected to the above constraint. Clearly, at 
values of y where the constraints are satisfied, 7 gives an 
upper bound on /obj(y)- The auxiliary polynomials Vjiy) 
and SOS /ij(y) then serve as the Lagrange multipliers in 
the relaxed optimization problem. 
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For a fixed degree of the above expression, this re- 
laxed optimization problem can be cast as an SDP in 
the form of Eq. ((9|) [28(. For the lowest order relaxation, 
the auxiliary polynomials Vj(y) and SOS /tij(x) are cho- 
sen such that degree of the expression in Eq. (|B2[) is no 
larger than the maximum degree of the set of polynomials 

/obj(y), / cq ,i(y), ■ • • , /cq,JV oq (y), /incq,l(y)> ■ • • , /incq,JV inoq (y); 

for a QCQP with no inequality constraints, this amounts 
to setting all the fii(y) to zero and all the Vj(y) to num- 
bers. 

For higher order relaxation, we increase the degree of 
the expression in Eq. (|B2|) by increasing the degree of the 
auxiliary polynomials. At the expense of involving more 
computational resources, a tighter upper of / bj(y) can 
then be obtained by solving the corresponding SDPs. 



APPENDIX C: LOWEST ORDER RELAXATION 
WITH OBSERVABLES OF FIXED TRACE 

To obtain a tighter upper bound by using the lowest 
order relaxation to Q18p . we found it most convenient to 
express (|18l) in terms of the real optimization variables, 



where X mn are defined in Eq. ([28 







-b®R' 



TtiA-\-n~iE 



2 \ -~b T ® R' T 



M m , 



1 = vec(L), r 



(C5a) 



y mn = tr (O m a n ) , n = 0, 1, . 



,d 2 - 1, 



(CI) 



which are just the expansion coefficients of each O m in 
terms of a set of Hcrmitian basis operators satisfying 
Eq. (jAlj) . The constraints can then be taken care of by 
setting each y m Q = z m j\fd. It is also expedient to ex- 
press the density matrix p in terms of the same basis of 
Hcrmitian operators 



d A -i 



([rA}i<Ji ® (Jo + [rs]icro <S> <Ji 



-*-d 2 

1=1 

d 2 -i 



where 



[R] i:j = tr(p cr, ®<7j) , (C2a) 
[r A ], = tr(p <Tj ® cr ), [r B ]j = tr(p a ® cTj); (C2b) 

r^, are simply the coherence vectors that have been 
studied in the literature [5^ |. 

We can thus incorporate the constraints (f2"5)) by ex- 
pressing the Lagrangian (|22|) as a function of the reduced 
set of variables 

y' T = [yu yi2 ■ ■■ yi^—i yn ■ ■■y mA +m B d?-i], ( C3 ) 

while all the y m o = z m j\fd are treated as fixed parame- 
ters of the problem. With this change in basis, and after 
some algebra, the Lagrangian can be rewritten as 



mA+rriB 



m— 1 

l 

7^ 



? 2 \ 



(i-r)V-y' T oy, 



k.l 



~J 2 

(C4) 



^ t B ® r B J 
and for i, j = 1, 2, . . . , d 2 — 1, 

[i]jm = 2z m A mj -, [t^Jfe = ^ b k iz [+mA , [tg]j = ^ &fc;^fc; 

i ft 

d 2 -l 

M m = ^ ^mnPn, [Pn]ij = ^ \ 
n=0 

(C5b) 

[<7i, <7j]+ = Oi<Jj + OjGi is the anti-commutator of <Ji and 

As before, we now maximize the Lagrangian (|C4[) 
over y' to obtain the corresponding Lagrange dual func- 
tion. The latter, however, is unbounded above unless 




V^(l-r) 20' 

for some finite t. The convex optimization problem dual 
to Eq. (|18p with fixed trace for each observables is thus 



mA+rriB 



minimize 



.2 \ 



a-.; 



(P 



( -2t -L (l T - r T ) \ 
subject to , ^ v ; > 0. (C7) 



APPENDIX D: NONLOCALITY CRITERION 



To derive the semianalytic criterion Eq. (|33[) , we now 

,2 i 

note that any choices of {A rnn }„ =0 that satisfy constraint 
(|C6[) will provide an upper bound on the corresponding 
Squ(p)- In particular, an upper bound can be obtained 
by setting 



6 n [Aa (8 m l + 6 m2 ) + \b (6m3 + <W)] (Dl) 



A, 



and solving for \a, that satisfy the constraint (|CJ6[> . 
With this choice of the Lagrange multipliers, and for 
quantum states with vanishing coherence vectors, the 
constraint (IC6I) becomes 



-2t 









-b®R' 
—b ® R T ^1 2 ®1«p_i J 



2ftl 2 ®l#-i 



> 0. (D2) 



where b and i?' are defined, respectively, in Eq. (|3Tj) and 
Eq. (|C2a[) . This, in turn is equivalent to 



-t > 0, 
-b®R 



-b®R' T 3Aj.l a ®l«p_i 



> 0. 



(D3a) 
(D3b) 
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Using Schur's complement [42|], the constraint (|D3b|) can 
be explicitly solved to give 



where s± is the largest singular value of the matrix R' . 
Substituting this and Eq. (|D3a[) into Eq. (IC7|) . and after 
some algebra, we see that (£>chsh)p for a quantum state p 
with vanishing coherence vectors cannot be greater than 



max 



Z 1; Z 2, Z Z, Z A 



\ 



n 



2d 2 



2d 2 



b k i- 



k,l 



Z kZl+2 
d 2 ' 



For p to violate the Bell-CHSH inequality, we must have 
this upper bound greater than the classical threshold 
value of 2 [c.f. Eq. Hence a sufficient condition for p 
to satisfy the Bell-CHSH inequality is given by Eq. 
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